Role of the Coulomb and the vector-isovector p potentials in the isospin asymmetry of 

nuclear pseudospin 



m 

O 

o 

(N 



R. Lisboa and M. Malheiro 

Instituto de Fisica, Universidade Federal Fluminense, 24210-340 Niteroi, Brazil 

P. Alberto 

Departamento de Fisica and Centra de Fisica Computacional, 
Universtdade de Coimbra, P-3004-516 Coimbra, Portugal 
(Dated: February 8, 2008) 

We investigate the role of the Coulomb and the vector-isovector p potentials in the asymmetry of 
the neutron and proton pseudospin splittings in nuclei. To this end, we solve the Dirac equation 
for the nucleons using central vector and scalar potentials with Woods-Saxon shape and Z and 
N — Z dependent Coulomb and p potentials added to the vector potential. We study the effect 
of these potentials on the energy splittings of proton and neutron pseudospin partners along a Sn 
isotopic chain. We use an energy decomposition proposed in a previous work to assess the effect 
of a pseudospin-orbit potential on those splittings. We conclude that the effect of the Coulomb 
potential is quite small and the p potential gives the main contribution to the observed isospin 
asymmetry of the pseudospin splittings. This isospin asymmetry results from a cancellation of the 
various energy terms and cannot be attributed only to the pseudospin-orbit term, confirming the 
dynamical character of this symmetry pointed out in previous works. 
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'. I. INTRODUCTION 
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i Pseudospin was introduced in the late 1960s dQl account for the quasidegeneracy of single- nucleon states with 
T-H ■ quantum numbers (n, £, j = £ + 1/2) and {n — 1, £ + 2, j = £ + 3/2) where n, £, and j are the radial, the orbital, and 
. the total angular momentum quantum numbers, respectively. These levels have the same "pseudo-orbital" angular 
' momentum quantum number, £ = £+1, and "pseudospin" quantum number, s = 1/2. 

, The pseudospin symmetry has been analyzed in nonrclativistic models by several works |^^|^. The subject was 
• revived in recent years when Ginocchio presented a relativistic interpretation for this symmetry showing that the 
I , quantum number £ is the orbital quantum number of the lower component of the Dirac spinor for spherical potentials. 
^ ■ Moreover, he showed that £ is a good quantum number in a relativistic theory for the nucleon with scalar S and vector 
^ ] V potentials with opposite signs and the same magnitude, i.e., S = S + V = 0. Actually, this condition can be relaxed 
^ ' to demand that only the derivative of E is zero 13], but, if E goes to zero at infinity, both conditions are equivalent. 
^ ' These findings are in agreement with the nuclear phenomenology obtained in relativistic mean field (RMF) theories, in 
• ^ view of the fact that the scalar and vector potentials in RMF for the nucleus 13, Tl'l , although large, cancel each 
^\ other to a great extent, giving a relatively small negative value for E on the nucleon mass scale .12.]. The structure 
of the radial nodes occurring in pseudospin levels as well as the nuclei wave function for them were also explained in 
recent works |3 0, ^| . Unfortunately, because in the nonrclativistic reduction of the Dirac equation E acts as the 
binding potential, it cannot be set to zero in nuclei, since in that case there would not be any bound state. 

So the question remains how the pseudospin symmetry is realized in nuclei, namely how some levels of pseudospin 
partners become degenerate or near degenerate. A partial answer to this question has been given in two previous works 
liallTl, in which were shown that, using Woods-Saxon potentials, this degeneracy not only depends on the value of 
the depth of the E potential but also on its characteristic extension (radius) and its surface steepness (diffusivity) . 
Furthermore, using an energy decomposition based on a Schrodinger-like equation for the radial lower component 
of the Dirac spinor, it was shown that the near-degeneracy of some pseudospin doublets arises from a cancellation 
between the terms in that decomposition, and that the pseudospin-orbit term is in general bigger than the splitting 
itself. This agreed with the findings of Refs. [isl IT^] in which the pseudospin-orbit coupling was shown to be a 
nonperturbative quantity. The conclusion was that pseudospin symmetry in nuclei has a dynamical character. 

The systematics found in Refs. flil ITtI for the pseudospin dependence on the shape of the nuclear mean fields 
allowed the explanation for the difference between the energy splittings of proton and neutron pseudospin partners. 
This was done by revealing the influence of the p potential Vp, which appears in RMF theories, on pseudospin-orbit 
interaction, favoring pseudospin symmetry for neutron levels as is found experimentally. However, the role of the 
Coulomb potential Vcoui for the proton splittings, which has the opposite sign of Vp, was not considered. A previous 
study of the isospin dependence of the pseudospin symmetry was done in Ref. pol | for the Zr and Sn isotopes, where 
it was shown that the pseudospin symmetry becomes better for exotic nuclei with a highly diffuse potential. 
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In this paper we examine in a more quantitative way this isospin asymmetry, by including exphcitly the Coulomb 
and the Vp pote ntials in the nuclear mean fields for the Sn isotopic chain and assess the role of each potential in that 
asymmetry [2lj . We will show that even though the pseudospin orbit term is very sensitive to the Coulomb potential 
Vcoui, the net effect of this potential in the proton pseudospin energy splitting is quite small. The main contribution to 
the asymmetry of the neutron and proton pseudospin splitting in nuclei comes from the vector-isovector Vp potential 
confirming the results found in Refs. |l6lll7|. 

The paper is organized as follows. In Sec. II we review briefly the formalism of the Dirac equation with scalar 
and vector potentials and the energy decomposition based on a Schrodinger-like equation for the lower component of 
the Dirac spinor. In Sec. Ill the neutron and proton mean field potentials in a Sn isotopic chain are parametrized 
by Woods-Saxon form as functions of A, N and Z including explicitly the Vp and Vcoui potentials. The variation of 
the pseudospin splittings with varying A and how Coulomb and p potentials affect those splittings, separating the 
contribution of each term in the energy decomposition, are presented in Sec. IV. In Sec. V we give a summary of the 
work and draw the conclusions. 



II. DIRAC EQUATION AND PSEUDOSPIN SYMMETRY 

The Dirac equation for a particle of mass m in an external scalar potential S and vector potential V is given by 

iJ* = [a ■ p + (3{m + S) + y]* = e * , (1) 

where a and (3 are the usual Dirac matrices. Defining A — V — S and the upper and lower components as = 
[(1 ± /3)/2]*, we have 

a-p*+ = (i; + 2TO-A)*_ (2) 

a-p*- = (3) 

where E = e — m is the binding energy. The Hamiltonian in Eq. ^ is invariant under SU(2) transformations when 
S = V or S = —V l2^ l23|. The second case corresponds to the realization of pseudospin symmetry. As referred 
before, this symmetry is related to the orbital quantum number of the lower component of the Dirac spinor This 
can be seen if one writes the corresponding decoupled second-order equations of Eqs. (0) and (O, assuming that the 
potentials S and A are radial functions: 

pH+ = - p^f ( -^ + -a-L)^+ + iE + 2m-A){E-j:)^+ (4) 
E + 2m — A \ or r J 

S' / d . 1 



p'^- = --^-^(^-^ + -(T-ij*_ + (S + 2m-A)(S-S)*_, (5) 

where the primes denote derivatives with respect to r and cr is the usual diagonal 4x4 spin matrix. The spinors 
can be separated into radial and angular parts, \E'_|_ = i G'k(?') . (6*, and = — ^^(r) ^^^^ . (6', </>), where k is 
the quantum number related to spin-orbit coupling and rrij is the eigenvalue of the Jz operator. Doing this separation 
in Eqs. and (O we arrive at the radial equations 



1 d f ^dGA + , A' / , , 1 



r 



dr \ dr J E + 2m — A 



K 



r 



G« + V + \+{E + 2m- A}(E - E)G„ = (6) 



1. /^,i£,VM±i)f. + ^fe + i^^ ,+,^5 + 2™-A)(i5-E)i..=0, (7) 



since one has cr ■ L^^^^ — —(1 ± k)^^^^ [2^. The term with 1 — k in Eq. {T)) is the pseudospin-orbit term. From 

that equation one sees that, should it be possible to set E' — 0, ^ would be a good quantum number. Since the sign 
of K determines whether one has parallel or antiparallel spin 

and £ = £ — k/\k\, one sees that pairs of states with k = —{£+1) and k = £ + 2 have the same £ = £+1, the quantum 
numbers of the pseudospin partners. 
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The second-order differential equation for Vt- [Eq. |SJ] can be written as a Schrodinger-type equation 



2m* 



1 



2 m,* E-J: 



d_ 

dr 



1 



El'. 



(9) 



where m* = (iJ + 2m — A)/2 is an energy- and r-dependent effective mass. This equation allows an energy decom- 
position to be performed by taking the expectation value of each term relative to the spinor 



where 
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In these expressions, P stands for the principal value of the integral. In Ref. ^3 it was shown that the pseudospin- 
orbit and Darwin terms ((Vpso) and (Vd), respectively) obey a sum rule that provides a numerical check for the 
computation of the principal values. The decomposition (|10|l will allow us to evaluate the importance of each term, 
and in particular of (Vpso)5 for the pseudospin splitting calculations shown in Sec. IV. 



III. PARAMETRIZATION OF THE NUCLEAR POTENTIAL FOR THE SN ISOTOPIC CHAIN 

In RMF theories, the main contributions to the isospin breaking in nuclei are given by the Lorcntz vector electrostatic 
(Coulomb) potential and the isovector- vector potential Vp coming from the interaction of p mesons with the nucleons 
|26|. These potentials modify the vector potential to 

Y ^ I Vc + Vp + Vboui for protons 

1 Vc — Vp for neutrons, ^ ' 

where Vc stands for the central (isoscalar) vector potential. Using the definition of E and A potentials we see that 
both are affected in the same way; i.e., 

E or A = / '-^'^ or A^) + Vp + Vboui for protons 

1^ (Ec or Ac) — Vp for neutrons, ^ ' 

where Yc — S + Vc and Ac = 14 — 5. As in Ref. 0,0], we use a Woods-Saxon form 

^^'^ = l + eJ[l-R)/a] ^''^ 
for Ec and Ac potentials with depth Uq, radius R, and diffusivity a to be determined later for each potential. 
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For Vp, we again use a Woods-Saxon shape, 

"'<'•' - l + .x,J-i;,)/a,l ^ 

for Vcoui we take the proton electrostatic potential energy in a uniform spherical charge distribution of radius Rc, 

r < Rc 
r > Rc- 



Vcoui(r)=<; ^\'"f^^^^ (19) 



47reo r 



In order to determine the effect of Vcoui a-nd Vp on pseudospin splittings, we solve the Dirac equation with the E and 
A potentials for a Sn isotopic chain, from A — 100 to ^4 = 170. Since there are almost no experimental data available 
for Sn neutron and proton single-particle energies, we fit the parameters of the central potentials, Vp [Eq. (|18|l ] and 
Vcoui [Eq. H19() ] to the S and A potentials given in Ref. [23| for each isotope. In that work a relativistic continuum 
Hartree-Bogoliubov calculation (RCHB) with a pairing interaction, which is an extension of the RMF theory, was 
performed to study tin isotopes. As a first step to fit the neutron potential I]„, we used a Woods-Saxon shape and 
found the corresponding parameters for each isotope. In the case of the proton potential Ep, we subtracted the 
Coulomb potential and fitted the nuclear part also by a Woods-Saxon form. In order to do that, we used in Eq. I|19() 
Rc = 1.20^1/3, which has the same A dependence of the nuclear radius. In Fig.nwe present the S„ and Ep potentials 
for each isotope obtained by our fits, together with the values of the corresponding Woods-Saxon parameters. These 
curves are very similar to the ones obtained by the RCHB calculation ^27\ which are also presented for comparison in 
Fig. n In the second step, using Eq. ifTCI) . we extracted Ec and Vp from those fits: 



(Ep - "t^coui + S„) (20) 



Vp = -(Ep-Vboui-E„) . (21) 

Finally, for each Sn isotope, we made a fit of Ec and Vp matching to Woods-Saxon shape. We found that the 
parameters correlate well either with mass number A or with proton and neutron number difference Z — N . This 
allowed us to have a general parametrization for the E potential for the whole isotopic chain: 

= TT rf°' p w 1 ± T- 7°' p w 1 + ^coui(r) (22) 

1 + cxp[(r - Rc)/ac\ 1 + exp[(r - Rp)/ap\ 

where the last term Vcoui exists only for protons. The parameters for Ec and Vp potentials in Eqs. (|20|l and H21|l . as 
functions of A, N, and Z, are 



Eq C — 


-69.94 MeV, 




(23) 


Rc = 


1.21^1/3 fm. 




(24) 




0.13^1/^ fm. 




(25) 


Vop = 


-[0.12(iV- Z) 


+ 3.87] MeV, 


(26) 


Rp — 


[om{N -z) + 


5.05] fm. 


(27) 


ap = 


[0.007(7V- Z) - 


f 0.27] fm. 


(28) 



The parametrization for Rc and Oc has a natural justification in view of the known A^^^ dependence of the nuclear 
radius. The Vop, Rp, and ap dependencies on iV — Z are also justified since they are proportional to the difference 
between proton and neutron densities. In the RCHB calculation the nuclear potentials for protons and neutrons for 
the case N = Z (^""Sn) are not exactly the same. The Coulomb potential changes the proton energy levels and 
because of the self-consistent RCHB calculation, the neutron energy levels are also changed in such a way that there 
is a net Vp potential. This effect was considered in the parametrization when we chose a value for Vop, Rp, and ap in 
the case oi N ^ Z a.s seen in Eqs. 

We show in Fig. |2the different terms that contribute to Ep and E„, respectively, for A = 150. In order to see the 
quality of our fit, we also present the respective E potentials obtained in the RCHB calculation [23]. From Fig. |2Ia), 
for the proton potential, we see that Vbp is around —10 MeV, while Vcoui is close to 20 MeV at the origin, so that the 
net effect at the origin increases Ec by ~ 10 MeV. In the case of the neutron potential in Fig.|2Ib), the net increase 
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1: (a) Neutron RCHB potentials for the Sn isotopic chain from Ref. [2^ . (b) Woods-Saxon fit for the neutron potentials 

bThe parameters R and a are in fermis and Eo in MeV.(c) Proton RCHB potentials for the Sn isotopic chain from 
. (d) Woods-Saxon fit for the proton potentials in (c). The units of the parameters are as in (b). 



at the origin turns out to be the same, because the only contribution is from Vp with a reversed sign compared to the 
proton case. Taking into account the dependence of the pseudospin splittings on the depth of the S potential only, 
one would expect that the proton and neutron splittings would be almost the same. However, as we will show later, 
this is not the case. The mechanism that generates the isospin asymmetry of pseudospin is more complex. 
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IV. THE ROLE OF Vcoui AND Vp IN THE ISOSPIN ASYMMETRY OF THE PSEUDOSPIN 

In this section, we single out the contribution of the Coulomb and vector-isovector Vp potentials to the pseudospin 
energy splitting for the Sn isotopic chain. These two potentials, as we already referred, are the main source of the 
isospin asymmetry of the nuclear pseudospin: the quasidegeneracy is better seen for neutrons than for protons. 

In order to identify the origin of this asymmetry and also to see if it has a direct relation to the pseudospin orbit 
potential, we used the energy decomposition given in Eq. The proton and neutron potentials obtained for the 

Sn isotopes with Woods-Saxon fit do not reproduce the small variations in the inner part of the RCHB potentials as 
we can see in Fig. Therefore, the deeper energy levels are not so well reproduced as the less bound ones. This 
should not be a problem insofar as pseudospin symmetry in nuclei is concerned, since it is seen experimentally that 
this symmetry is in general much better satisfied for the levels near the Fermi sea. To perform our analysis we chose 
the [2dc^i2 — I37/2] pseudospin partner, which is close to the Fermi sea. If we had selected any other doublet near the 
Fermi sea, the conclusion of the foregoing analysis would have been the same. 

The contribution of each term of Eq. l(Tn|l to the pseudospin energy splitting for the [2^5/2 — 137/2] doublet, with 
and without the inclusion of the Coulomb and vector-isovector Vp potentials, is shown in Fig. O for the Sn isotopes 
from A — 130 to A ^ 150. The Sn isotopes with mass number in the range A =100-130 were not considered because 
we found that in some cases the above levels were not bound when we excluded the Vp potential. The effect of Vcoui 
or Vp on the energy splitting of the pseudospin partners is presented in Fig. 0] Let us first analyze the effect of the 
Coulomb potential. In Fig. |3[c) we see that the contribution from the pseudospin orbit term {Vps) to the proton 
energy splitting changes considerably when we exclude the Coulomb potential. This contribution becomes smaller 
(— 1 MeV) and essentially the same as the one for the neutron energy splitting. Thus, if we link the pseudospin 
symmetry only to the pseudospin orbit term and do not take in account the contribution of the other terms, we will 
be forced to conclude that the Coulomb potential is the essential source of the isospin asymmetry seen in the nuclear 
pseudospin. However, as it has already been pointed out in a previous work ^3|> the pseudospin degeneracy comes 
from a significant cancellation among all the terms in Eq. (jlOII and not only by the (Vps) interaction, as can be seen 
in Figs. Ola) and|3Jd). In fact, the contribution for the proton energy difference from kinetic {p^ /2m*) and potential 
(S) energies increases approximately the same for each of these terms (~ 0.5 MeV) when we exclude the Coulomb 
interaction. These effects almost cancel the changes produced in the pseudospin orbit term (Vps). Thus, the small 
net effect produced by the exclusion of the Vcoui in the proton pseudospin energy splitting presented in Fig. 0] is 
essentially due to the changes in the Darwin term, which are quite small as shown in Fig. |3^b). The fact that the 
Coulomb potential docs not explain the isospin asymmetry reflects the dynamic nature of the nuclear pseudospin 
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symmetry. It comes from the contribution of all the terms to the energy splitting and cannot be attributed only to 
the pseudospin-orbit term {Vps), which in this case is strongly affected by (Vcoui)- 
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FIG. 3: The effect of Vcoui and Vp in each term of Ea. ljlO^ that contributes to the neutron and proton pseudospin energy 
splittings for the [2^5/2 — I57/2] doublet in the Sn as a function of the mass number A. (a) Kinetic term, (b) Darwin term, (c) 
pseudospin-orbit term, (d) potential E term. 



Let us now analyze the effect of the vector-isovector potential Vp in this isospin asymmetry. First, we see from 
Fig-EIc) that the contribution of the term (Vps) to the neutron and proton energy difference remains almost unchanged 
when we exclude the Vp potential. The asymmetry between protons and neutrons, as we see in Figs.|3Ia) and|3|[d), 
comes essentially from the effect of Vp in the kinetic and potential terms. When we exclude this potential in the case 
of the neutrons, the contribution of {p^ /2m*) and (S) terms are bigger (the E potential is deeper), thus increasing the 
neutron pseudospin energy difference as we can see in Fig. ^ In the case of the protons we see an opposite behavior 
for these two terms, since the depth of the E potential is now smaller {Vp is attractive for protons), decreasing the 
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FIG. 4: The effect of Vcoui and Vp in the neutron and proton pseudospin energy 
sphttings for the [2^5/2 — I57/2] doublet in the Sn as a function of A. 



proton pseudospin energy difference. Thus, as we see in Fig. ^ the neutron and proton pseudospin energy sphttings 
approach each other becoming almost degenerate. This analysis allow us to conclude that the origin of the proton 
and neutron asymmetry comes mainly from the vector-isovector Vp potential. 

This conclusion is in agreement with previous works that have also investigated this isospin asymmetry |l6l Il7j | . 
Our analysis reveals that the effect of Vp in the pseudospin potential in Eq. H10() shows up mainly in an indirect way, 
i.e., through the change induced in the lower component of the wave function "if-. This effect has the result that, 
contrary to what could be expected, the contribution of the expectation value (Vps) of the potential to the pseudospin 
energy splitting, either for protons or neutrons, is almost insensitive to the vector-isovector potential. This fact again 
manifests the dynamical nature of the nuclear pseudospin symmetry. The isospin asymmetry analyzed here is then 
a consequence of the effects of Vp potential in all the terms that contribute to the proton and neutron pseudospin 
energy splittings. 



V. CONCLUSIONS 



We have investigated in a quantitative way the effect of the Coulomb and the vector-isovector p potentials in the 
proton and neutron asymmetry seen in the nuclear pseudospin. To do this analysis we performed a mean-field model 
calculation with Woods-Saxon potentials fitted to the proton and neutron potentials obtained by a sophisticated 
self-consistent RCHB calculation for the Sn isotopes . We found that the Woods-Saxon parameters correlate well 
with the mass number A and the proton and neutron number difference N — Z, which allowed us to have a general 
parametrization of the binding potential for the whole isotopic chain. In this general parameterization of the potential 
we have separated explicitly the Vcoui and Vp potentials. These two potentials are the main source of the isospin 
asymmetry of the nuclear pseudospin. 

In order to identify the origin of this asymmetry we analyzed the effect of those potentials separately in the proton 
and neutron pseudospin energy splitting along the Sn isotopic chain. This analysis was done taking in account the 
effect of these potentials in each of the terms that contribute to the pseudospin energy splitting. This decomposition 
of the energy helped us to uncover the mechanism behind this isospin asymmetry. We showed that the isospin 
asymmetry seen in the pseudospin energy splitting cannot be attributed to the pseudospin orbit term. It comes from 
a strong cancellation among all the terms that contribute to the energy in agreement with the dynamical nature of 
these symmetry pointed out before in 0, 0, 0, 0| . 
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We explicitly showed that even though the pseudospin orbit term (Vps) is strongly affected by the Coulomb potential 
Vcouh the net effect of this potential on the difference of proton and neutron pseudospin energy is quite small because 
of the effects of Vcoui in all the other terms. 

The main conclusion of our paper is that the vector-isovector Vp potential gives the main contribution to the 
observed isospin asymmetry of the nuclear pseudospin, as stated in previous works [l^ Il7|. However, it has little 
effect on the contribution of the pseudospin orbit interaction term < Vps > to the pseudospin energy splitting, which 
means that the changes of the pseudospin potential and of the wave function F^^ir) caused by Vp seem to cancel 
each other. Therefore, the role of Vp in the pseudospin isospin asymmetry comes from its effects in the other terms 
contributing to the single-particle energy, especially in the differences between the kinetic and potential energy terms 
of the pseudospin partners. This important result shows that the dynamical character of pseudospin symmetry and 
in particular the nonperturbative nature of the pseudospin potential lie at the heart of the mechanism responsible for 
the differences of proton and neutron pseudospin energy splittings observed in nature. 
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